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A  previously  published  computational  multi-phase  model  of  a  polymer-electrolyte  membrane  fuel  cell 
cathode  has  been  extended  in  order  to  account  for  the  anode  side  and  the  electrolyte  membrane.  The 
model  has  been  applied  to  study  the  water  balance  of  a  fuel  cell  during  operation  under  various  humidifi¬ 
cation  conditions.  It  was  found  that  the  specific  surface  area  of  the  electrolyte  in  the  catalyst  layers  close 
to  the  membrane  is  of  critical  importance  for  the  overall  water  balance.  Applying  a  high  specific  elec¬ 
trolyte  surface  area  close  to  the  membrane  (a  water-uptake  layer)  can  prevent  drying  out  of  the  anode  and 
flooding  at  the  cathode  while  the  average  membrane  water  content  is  only  weakly  affected.  The  results 
also  indicate  that  in  contrast  to  common  presumption  membrane  dehydration  may  occur  at  either  anode 
or  cathode  side,  entirely  depending  on  the  direction  of  the  net  water  transport  because  the  predominant 
transport  mechanism  is  diffusion.  Consequently,  operating  conditions  with  a  high  net  water  transport 
from  anode  to  cathode  should  be  avoided  as  it  is  important  to  keep  the  cathode  catalyst  layer  well  humid¬ 
ified  in  order  to  prevent  high  protonic  losses.  Addition  of  the  micro-porous  layer  did  not  affect  the  overall 
water  balance  or  membrane  water  content  in  our  study. 

©  201 1  Elsevier  B.V.  Ah  rights  reserved. 


1.  Introduction 

A  fundamental  understanding  of  the  water  balance  of  a  fuel  cell 
during  operation  is  crucial  for  improving  the  cell  performance  and 
durability.  It  is  well  known  that  an  improper  water  management 
will  lead  to  excessive  flooding  of  the  cathode  and  presumably  to 
drying  out  of  the  anode  as  water  is  being  dragged  from  anode 
to  cathode  along  with  the  protonic  flux  through  the  membrane, 
a  process  commonly  named  electro-osmotic  drag  (EOD).  The  EOD 
coefficient  has  been  experimentally  determined  to  be  around  unity 
by  Zawodzinski  et  al.  [1]  and  Ye  and  Wang  [2],  meaning  that  one 
water  molecule  per  proton  is  dragged  from  anode  to  cathode.  An 
EOD  of  unity  already  means  that  twice  the  amount  of  product  water 
created  at  the  cathode  is  dragged  from  anode  to  cathode,  and  the 
conventional  means  of  providing  sufficient  water  to  the  anode  to 
prevent  drying  out  is  to  humidify  the  incoming  anode  gas  stream 
[3]  and  to  allow  for  back  diffusion  of  water  from  cathode  to  anode. 
To  this  end  the  diffusivity  of  water  through  the  membrane  has 
been  measured  by  several  research  groups  (e.g.  [4]).  Recently,  Mon¬ 
roe  et  al.  [5]  have  conducted  experiments  on  a  plane  membrane 
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and  found  that  for  thin  (<50  pun)  membranes  as  they  are  being 
used  nowadays  the  overall  water  transport  may  not  be  controlled 
by  diffusion  inside  the  membrane  phase,  but  also  by  the  absorp¬ 
tion/desorption  kinetics  of  water  to/from  the  electrolyte  phase. 
Ge  et  al.  [6]  have  in  turn  measured  that  water  absorption  to  the 
membrane  is  around  one  order  of  magnitude  slower  than  desorp¬ 
tion,  and  their  experiments  were  conducted  on  the  operating  fuel 
cell. 

In  this  work  we  have  extended  a  previously  published  model 
[7,8]  in  order  to  account  for  the  anode  side  and  the  polymer- 
electrolyte  membrane.  Water  transport  across  the  membrane  is 
modeled  using  a  constant  EOD  of  unity  and  a  diffusion  coefficient  of 
water  in  the  electrolyte  phase  taken  from  the  literature.  The  absorp¬ 
tion  coefficient  was  also  taken  from  the  literature,  but  in  our  study 
the  specific  surface  area  of  the  electrolyte  in  the  catalyst  layer  to 
uptake  the  water  was  varied  over  one  order  of  magnitude.  It  will  be 
shown  that  this  parameter  is  very  important  for  the  overall  water 
balance.  We  also  confirm  that  for  thin  membranes  the  diffusion  of 
water  in  the  electrolyte  may  not  be  rate  determining. 

For  the  water  balance  it  is  desirable  that  the  fuel  cell  cathode 
does  not  become  flooded  and  the  anode  does  not  become  dry.  In 
addition,  the  water  content  X  of  the  electrolyte  membrane  should 
be  as  high  as  possible  in  order  to  reduce  ohmic  losses.  A  layer  with 
a  high  specific  surface  area  a  of  the  electrolyte  to  facilitate  suffi¬ 
ciently  fast  water  absorption  is  called  a  water-uptake  layer  in  our 
work,  and  it  is  found  that  such  a  layer  may  substantially  increase 
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Nomenclature 

a  water  activity 

C  mass  concentration  (kg  m-3 ) 

D  diffusion  coefficient  (m2  s-1 ) 

EW  equivalent  weight  of  dry  membrane  (=1 .1  kg  mol-1 ) 

F  Faraday’s  constant  (=96,485  C  mol-1 ) 

fv  specific  volume  fraction  of  water  in  the  membrane 

fi  specific  volume  fraction  of  water  in  the  liquid- 

equilibrated  membrane 
i"  current  density  (A  cm-2 ) 

I<  permeability  (m2) 

ka,  kd  kinetic  absorption/desorption  coefficient  (m  s-1 ) 
krd  relative  permeability 

M  molecular  weight  (kg  mol-1 ) 

nd  drag  coefficient 

P  pressure  (Pa) 

5  source  term 

s  saturation 

rdrag  effective  (net)  drag  coefficient 

T  temperature  (I<) 

U  velocity  vector  U  =  (u,  v,  w) 

Vm  partial  molar  volume  of  dry  membrane 
Vw  partial  molar  volume  of  membrane  water 

Y  dimensionless  mass  fraction 

Greek  symbols 

a  specific  electrolyte  surface  area  (m2  m-3 ) 

s  porosity 

0  contact  angle 

A  membrane  water  content 

v  molecular  viscosity  (m2  s-1 ) 

fi  dynamic  viscosity  (kg  m- 1  s-1 ) 

p  density  (kg  m-3) 

pmem  ^ry  membrane  density  (=2000  kg  m-3 ) 
r  product  of  density  and  diffusion  coefficient 

Subscripts 

a  absorption 

d  desorption 

equil  equilibrium  value 

w  (membrane)  water 

Superscripts 
mem  membrane 

Abbreviations 

BP  bipolar  plate 

CL  catalyst  layer 

CFD  computational  fluid  dynamics 

EOD  electro-osmotic  drag 

GDL  gas  diffusion  layer 

MEA  membrane  electrode  assembly 

MPL  micro-porous  layer 

PEMFC  polymer  electrolyte  membrane  fuel  cell 

RH  relative  humidity 


the  limiting  current  density  when  operating  the  cell  under  dry  inlet 
conditions.  In  practice,  such  a  layer  can  be  applied  by  increasing 
the  mass  fraction  of  electrolyte  phase  inside  the  CL  in  the  vicin¬ 
ity  of  the  membrane  as  an  integral  part  of  a  multi-layer  electrode 
[9]. 


2.  Literature  review 

There  are  a  number  of  publications  on  water  balance  in  a  fuel 
cell  employing  both  experimental  and  numerical  methods.  An 
overview  of  water  balance  issues  and  experimental  efforts  was 
given  by  Dai  et  al.  [10],  and  Baschuk  and  Li  gave  a  summary  of 
ion  and  water  transport  phenomena  in  membranes  [11]. 

Theoretical  studies  employed  either  analytical  models  (e.g. 
[4,12,13])  or  the  methods  of  computational  fluid  dynamics  (CFD) 
(e.g.  [14]).  While  publications  of  CFD  models  for  PEMFC  abound 
in  the  literature  [15]  there  are  only  few  detailed  studies  pub¬ 
lished  on  the  water  balance  during  operation.  This  is  particularly 
unfortunate  because  water  balance  is  one  of  the  few  non-intrusive 
experiments  that  can  be  reliably  conducted  and  thus  used  for 
model  validation.  Luo  et  al.  [16]  presented  probably  the  most 
detailed  fuel  cell  CFD  model  to  date,  but  did  not  publish  a  detailed 
study  on  water  balance  in  a  PEMFC.  The  study  published  by  Lee 
et  al.  [14]  was  based  on  the  commercial  fuel  cell  module  by  CD 
Adapco,  which  is  essentially  single  phase  and  contains  only  a  sim¬ 
plified  membrane  model.  Both  models  do  not  account  for  the 
sorption/desorption  kinetics  of  water  into  and  out  of  the  electrolyte 
phase  but  instead  assume  equilibrium  between  the  electrolyte 
phase  and  the  adjacent  gas  phase  in  terms  of  water  content. 
Wu  et  al.  [17]  and  Ye  and  Nguyen  [18]  published  CFD  model¬ 
ing  studies  that  accounted  for  the  sorption/desorption  kinetics 
but  assumed  constant  proportionality  values  and  did  not  inves¬ 
tigate  the  effect  of  the  kinetics  and  specific  surface  area  of  the 
electrolyte.  Liu  et  al.  [19]  conducted  experiments  using  a  30p,m 
Gore  membrane  and  obtained  a  distributed  net  drag  coefficient  via 
postprocessing.  The  calculated  net  drag  was  low  and  even  nega¬ 
tive  due  to  strong  back  diffusion.  Finally,  Wu  et  al.  [17]  published 
a  three-dimensional  multi-phase  CFD  model  that  accounted  for 
membrane  water  transport  and  included  sorption  and  desorption 
terms.  Unfortunately,  they  did  not  publish  a  detailed  water  balance 
study  either. 

Experimental  work  on  water  transport  across  the  membrane 
was  reported  by  Janssen  and  Overvelde  [20],  Cai  et  al.  [21]  and 
Karan  et  al.  [22],  among  others.  Janssen  and  Overvelde  applied 
various  humidification  conditions  at  both  anode  and  cathode  side 
corresponding  to  either  fully  saturated  or  entirely  dry  gases.  The 
resulting  effective  drag  coefficient  was  positive  (water  transport 
from  anode  to  cathode)  for  most  cases.  However,  experiments  con¬ 
ducted  on  a  Nation  122  membrane  (50|jim  thickness)  also  show 
that  the  net  drag  can  be  negative  for  such  thin  membranes  under 
fully  humidified  conditions  [20]. 

Karan  et  al.  [22]  performed  experiments  with  an  (inlet-)  RH 
gradient  between  anode  and  cathode  and  found  that  the  drag  coef¬ 
ficient  was  positive  for  the  cases  investigated  (A/C  100%/60%  and 
60%/100%  with  and  without  MPL,  respectively).  Yan  et  al.  [23]  took 
water  balance  data  at  various  operating  conditions  and  found  that 
depending  on  the  RH  of  the  feed  gases  the  net  drag  could  be  either 
positive  or  negative.  Ye  and  Wang  [24]  measured  a  positive  drag 
coefficient  over  a  wide  range  of  operating  conditions.  Overall  it  can 
be  stated  that  experiments  were  published  only  for  a  few,  chosen 
operating  and  material  conditions. 

It  is  important  to  realize  that  most  of  above  studies  were  con¬ 
ducted  using  relatively  low  stoichiometric  flow  ratios  of  2-3.  While 
this  provides  a  better  match  with  realistic  fuel  cell  operating  con¬ 
ditions  it  is  bound  to  yield  three-dimensional  effects  in  any  given 
set-up.  Consequently  the  local  net  drag  will  vary  from  inlet  to  outlet 
and  depend  on  co-  and  counter  flow  conditions.  The  measured  lump 
value  for  the  net  drag  is  an  integral  value  of  the  local  values,  which 
makes  these  data  difficult  to  compare.  While  above  experiments 
provided  information  about  net  water  transport  under  realistic 
operating  conditions  they  are  not  very  suitable  for  improving  the 
fundamental  understanding. 
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Domain  V  Porous  membrane  for 
liquid  water  transport 


Domain  IV 


MEA  -  Electrolyte 
phase  for  dissolved 
water  diffusion  and 
EOD 


Bipolar  plate  - 

Domain  III  Heat  conducting 
solid 


Gas  and  liquid  phase 
phase  domain  with 
Domain  I  impermeable 

membrane  as  heat 
conducting  solid 


Domain  II 


Bipolar  plate  - 
Heat  conducting 
solid 


Darcy’s  law  for  liquid  water  permeation 
through  the  porous  membrane 
Boundary  conditions:  liquid  phase 
pressure  in  cells  from  Domain  I 

Diffusion  of  dissolved  water  in 
electrolyte  phase 

Source  and  sink  terms  in  exchange 
with  Domain  I  due  to  absorption/ 
desorption  in  CL 

Heat  conduction 

Specified  cell  temperature  at  outer  wall 
Temperature  boundary  condition  equal 
to  adjacent  cell  temperature  in 
Domain  I 

Source-  sink  terms  due  to  thermal 
contact  resistance  at  GDL  interface 

Navier-Stokes  equations  for  gas  phase 
and  liquid  phase 

Darcy  s  law  with  relative  permeabilities 
in  porous  media 
Phase  change  of  water 
Leverett  function  to  couple  capillary 
pressure  and  saturation 
Sink  and  source  terms  in  CL  due  to 
electrochemical  reaction 
Temperature  field  in  heat  conducting 
membrane 

Heat  conduction 

Specified  cell  temperature  al  outer  wall 
Temperature  boundary  condition  equal 
to  adjacent  cell  temperature  in 
Domain  I 

Source-  sink  terms  due  to  thermal 
contact  resistance  at  GDL  interface 


Fig.  1.  Two-dimensional  cut  through  the  computational  domains  and  summary  of 
physical  models  applied  [46]. 


Ideally,  water  balance  measurements  should  be  conducted  at 
a  high  stoichiometric  flow  ratio  in  order  to  limit  the  problem  to 
be  two-dimensional  (no  down-the-channel  variation).  If  in  addi¬ 
tion  the  channel-rib  spacing  is  chosen  very  fine  we  are  left  with 
a  quasi  one-dimensional  problem  which  will  proof  to  be  complex 
enough  for  water  balance  studies.  Experiments  conducted  in  the 
three-dimensional  mode  are  of  limited  value  for  model  validation 
purposes  because  of  the  large  number  of  material  parameters  and 
uncertainties  involved. 

While  detailed  experiments  on  water  balance  can  be  very  costly, 
a  simulation  with  this  model  took  approximately  2h  per  data 
point  on  a  laptop  computer  with  2.2  GHz,  which  allowed  for  a  very 
detailed  case  study.  The  drawback  is  that  no  CFD  model  of  a  fuel  cell 
has  been  extensively  validated.  Models  are  currently  being  used  to 
shed  additional  insight  into  experiments  conducted  and  can  also  be 
used  to  suggest  new  experiments  that  verify  some  of  the  modeling 
findings. 

3.  Model  description 

The  model  is  based  on  the  formerly  commercial  CFD  package 
CFX-4.4  ( ANSYS  Inc.).  The  advantage  of  this  code  is  the  implementa¬ 
tion  of  a  multi-fluid  model  and  the  possibility  to  customize  the  code 
via  an  extensive  suite  of  FORTRAN  user  subroutines.  In  this  section 
we  will  give  a  description  of  the  physical  model.  Fig.  1  shows  a  two- 
dimensionaly-z-cut  of  our  computational  domains.  The  x-direction 
is  defined  as  the  down-the-channel  direction. 

It  can  be  seen  that  we  have  utilized  five  separate  computa¬ 
tional  domains  to  model  the  various  transport  processes  that  occur 
in  a  fuel  cell.  These  domains  communicate  with  each  other  via 
implemented  source  and  sink  terms  as  will  be  described  below. 
It  is  important  to  realize  that  Domain  IV  is  an  exact  copy  of  the 
membrane-electrode  assembly  (MEA)  cells  from  Domain  I  so  that 
there  is  a  direct  cell-to-cell  correspondence.  Likewise  Domain  V , 
which  consists  of  the  membrane  only,  is  an  exact  copy  of  the 
membrane  from  Domain  /,  but  while  the  membrane  in  Domain  I 
is  represented  by  a  non-permeable  heat-conducting  solid,  Domain 
V  is  a  porous  region  where  the  Darcian  transport  of  liquid  water 
through  the  membrane  is  modeled  with  appropriate  pressure 
boundary  conditions  and  employing  a  local  hydraulic  permeabil¬ 


ity  that  is  a  function  of  the  membrane  water  content  calculated  in 
Domain  IV.  More  details  of  the  purposes  and  equations  applied  in 
the  different  computational  domains  are  given  below. 

3.1.  Domain  I 

This  is  the  main  computational  domain  where  we  solve  for  the 
gas  phase  and  liquid  phase  flow  in  the  channels  and  porous  media: 
the  gas  diffusion  layer  (GDL),  micro-porous  layer  (MPL)  -  if  included 
-  and  the  catalyst  layer  (CL).  For  each  thermodynamic  phase  a 
complete  set  of  Navier-Stokes  equations  (plus  energy  equation)  is 
solved.  In  addition  species  conservation  equations  for  oxygen  and 
water  vapor  at  the  cathode  side  and  hydrogen  and  water  vapor  at 
the  anode  side  are  solved.  Phase  change  of  water  is  accounted  for  at 
both  anode  and  cathode  side.  Boundary  conditions  are  prescribed 
at  the  channel  inlets  in  terms  of  gas  composition  and  temperature 
(no  liquid  enters  here),  while  we  specify  pressure  boundary  condi¬ 
tions  at  the  outlet.  In  the  porous  media  the  momentum  equations 
reduce  to  Darcy’s  law  with  a  different  relative  permeability  of  each 
phase.  The  liquid  phase  is  transported  here  by  capillary  forces  and 
the  capillary  pressure  versus  liquid  saturation  is  specified  using 
the  Leverett  equation  [25,26]  but  only  applied  to  the  reducible  sat¬ 
uration  and  thus  allowing  for  simulating  a  porous  medium  with 
a  mixed  wettability  (hydrophilic  and  hydrophobic  pores).  In  our 
model  a  different  irreducible  saturation  can  be  specified  for  every 
layer,  for  example  the  catalyst  layers  contain  a  much  larger  fraction 
of  hydrophilic  pores  than  the  MPL.  We  have  previously  argued  that 
the  shape  of  the  Leverett  function  should  be  closely  related  to  the 
integrated  pore-size  distribution  [7]  and  this  could  conceivable  be 
captured  by  adjusting  the  coefficients  in  the  cubical  Leverett  equa¬ 
tion.  Finally,  a  different  effective  liquid  water  contact  angle  0  can  be 
specified  for  every  porous  region.  This  contact  angle  applies  only  for 
the  hydrophobic,  water  transporting  pores.  Table  1  lists  important 
material  parameters  employed  in  the  current  study.  In  our  previous 
study  it  was  found  that  the  irreducible  saturation  predominantly 
dictates  the  saturation  level  encountered  in  every  layer.  Specifying 
different  material  properties  means  that  we  will  encounter  a  sat¬ 
uration  jump  between  every  adjacent  layer  because  the  capillary 
pressure  has  to  be  continuous  across  each  interface  [27]. 

Inside  the  three-dimensional  cathode  CL  there  is  oxygen  con¬ 
sumption  and  production  of  liquid  water  according  to  a  specified 
local  current  density  that  depends  on  the  local  oxygen  concen¬ 
tration.  A  detailed  electrochemical  model  is  not  yet  implemented. 
Along  with  the  water  production  there  is  heat  production  due  to  an 
assumed  local  overpotential  as  described  in  reference  [7].  To  sum¬ 
marize,  Domain  I  is  utilized  to  calculate  the  gas  phase  and  liquid 
phase  flow  fields  in  the  gas  flow  channels  and  porous  media.  Anode 
side  and  cathode  side  are  separated  by  an  impermeable  membrane 
that  only  allows  for  heat  conduction.  Species  equations  are  solved 
for  oxygen,  water  vapor  and  hydrogen  while  nitrogen  is  assumed 
as  background  fluid.  There  is  no  oxygen  or  nitrogen  at  the  anode 
and  no  hydrogen  at  the  cathode.  The  exchange  of  water  between 
anode  and  cathode  is  calculated  in  Domain  IV  and  Domain  V. 

3.2.  Domains  II  and  III 

These  domains  represent  the  bipolar  plates  (BPs)  and  currently 
they  are  only  used  to  calculate  heat  conduction.  Consequently  the 
only  equation  of  interest  here  is  the  energy  equation  to  yield  the 
temperature  distribution  inside  the  bipolar  plates.  Boundary  con¬ 
ditions  have  been  prescribed  at  the  outside  walls  in  form  of  the 
desired  cell  operating  temperature.  The  computational  cells  adja¬ 
cent  to  the  flow  channels  communicate  with  cells  in  Domain  I 
by  applying  the  temperature  in  the  adjacent  flow  channel  cell 
as  the  boundary  condition  for  Domains  II  and  III  and  vice  versa, 
which  makes  for  a  very  robust  numerical  scheme.  A  thermal  con- 
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Table  1 

Porous  material  properties  applied  in  the  current  study. 


Region  name 

Irreducible  saturation  [-] 

Contact  angle  [°] 

Permeability  [m2] 

Porosity  [-] 

Tortuosity  [-] 

GDL 

0.20 

120 

1.0E-12 

0.75 

4 

MPL 

0.10 

130 

0.1E-12 

0.75 

6 

CL 

0.30 

110 

0.1E-12 

0.75 

4 

tact  resistance  at  the  interface  between  the  BP  and  the  GDL  has 
been  implemented  to  account  for  contact  losses.  In  the  future  this 
domain  will  also  be  used  to  calculate  the  electrical  field  equation. 

3.3.  Domain  IV 


r\mem 

UW 


x 


10-10  exp  1 241 6  (^-L  -  j)  |  (2.563  -  0.33A 


+0.0264A2  -  0.000671  A3) 


(6) 


This  domain  represents  the  electrolyte  phase  of  the  MEA,  and  it 
accounts  for  the  transport  of  water  between  anode  and  cathode 
side.  The  only  equation  of  interest  solved  here  is  the  conserva¬ 
tion  of  dissolved  water  in  the  electrolyte  phase  [28].  Thus  a  pure 
diffusion  equation  is  solved  with  source  terms  inside  the  CL  due 
to  sorption/desorption  of  water.  The  electro-osmotic  drag  coeffi¬ 
cient  is  assumed  constant  in  this  work  which  leads  to  an  additional 
non-zero  source  term  for  the  electrolyte  water  inside  the  CL.  The 
diffusion  coefficient  is  itself  a  function  of  membrane  water  content 
and  temperature,  and  it  has  been  taken  from  the  literature  [4]. 

Mathematically  the  transport  equation  that  is  solved  here  is  the 
standard  convection-diffusion  equation  in  CFX-4  [28]: 

V  •  (yoUYw)  -  v  •  (rwVYw)  =  S  (1) 


By  comparison  to  Eq.  (5)  Luo  et  al.  [16]  solved  the  following 
equation  in  their  membrane  model: 

v  (gv„..v*)  Mw  -  V  ■  (4)  M„  +  V  (*=**)  .  o 

(7) 

The  second  term  in  their  equation  accounts  for  the  water  trans¬ 
port  due  to  electro-osmotic  drag.  It  can  be  easily  shown  that  for 
a  constant  drag  coefficient  this  term  is  zero  inside  the  membrane 
and  non-zero  inside  the  CL  [16,29].  In  our  model  this  is  accounted 
for  by  a  source  term  similar  to  Gurau  et  al.  [28]: 

“  W  V  '  =  S^sorp/desorp)  ~  V  •  (nd- J  Mw  (8) 


where  p  is  the  density,  U  is  the  bulk  velocity  vector,  rw  is  the 
product  of  density  and  the  diffusion  coefficient,  Yw  the  water  mass 
fraction  in  the  membrane  and  S  is  a  source  term.  The  convective 
term  (first  term)  inside  the  membrane  was  set  to  zero  as  there  is  no 
net  velocity  of  the  electrolyte  phase,  which  leaves  only  the  diffusion 
and  source  term.  The  mass  fraction  Yw  can  be  replaced  by  the  mass 
concentration  Cw  divided  by  the  density  of  the  membrane  pmem- 
This  yields  [16]: 
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Neglecting  the  swelling  of  the  membrane  the  density  terms  can 
be  cancelled  against  each  other.  The  correlation  between  the  con¬ 
centration  of  water  and  the  dimensionless  water  content  A  was 
given  by  Springer  et  al.  [4]: 
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While  the  dimensionless  mass  fraction  solved  for  by  CFX-4  is 
that  value  divided  by  the  membrane  density: 

=  0-016364  X  A  (4) 

This  yields  for  Eq.  (2): 


The  third  expression  in  Eq.  (7)  accounts  for  liquid  phase  trans¬ 
port  due  to  convection.  In  our  model  this  part  is  treated  in  Domain 
V  (see  below).  So  in  contrast  to  Luo  et  al.  we  treat  the  water  dis¬ 
solved  in  the  membrane  phase  differently  from  the  water  that  is 
permeating  through  the  membrane  by  pressure  forces. 

An  alternative  description  of  membrane  water  transport  was 
given  by  Janssen  [30]  who  published  a  model  for  water  and  pro¬ 
tonic  transport  through  the  membrane  using  concentrated  solution 
theory.  A  similar  approach  was  used  by  Weber  and  Newman  [31  ] 
and  Fimrite  et  al.  [32]. 

The  source  term  due  to  sorption/desorption  of  water  depends 
on  the  sorption  kinetics,  the  specific  surface  area  of  the  electrolyte 
inside  the  CL  and  a  driving  force  for  absorption/desorption  due  to 
deviation  from  the  equilibrium  value  of  kequu  (or  Cwequii )  [28]: 

^X(sorpldesorp)  =  OL  x  ka  X  (Cwequ //  —  Cw)  x  Mw  (9) 

In  case  of  water  absorption  to  the  electrolyte  phase  there  is 
a  source  term  for  the  water  in  the  electrolyte  phase  and  a  cor¬ 
responding  sink  term  for  the  gas  phase  water  in  Domain  I.  The 
specific  surface  area  of  the  electrolyte  phase  in  the  CL,  a  in  [m2  m-3  ], 
is  central  for  the  further  study.  Unfortunately,  no  detailed  data 
is  available  in  the  literature,  and  we  will  leave  this  parameter 
adjustable.  The  kinetic  sorption  coefficient  has  been  determined 
by  Ge  et  al.  [33]: 
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This  group  also  measured  the  desorption  rate  to  be  [33]: 


-Pmem^-(D™mV\)  =  S  (5) 

where  pmem  is  the  density  of  the  dry  membrane  (2000  kg  m-3)  and 
EW  its  equivalent  weight  (1100  kg  kmol-1 ).  A  very  similar  transport 
equation  is  used  by  Gurau  et  al.  [28].  It  must  be  considered  though 
that  CFX-4  can  only  solve  a  transport  equation  in  the  form  of  Eq.  (1 ) 
with  the  restriction  that  Yw  is  smaller  than  or  equal  to  unity.  For 
the  diffusivity  of  water  inside  the  membrane  Dw,mem  we  used  the 
expression  suggested  by  Springer  et  al.  [4]: 
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However,  in  our  steady-state  model  the  absorption  rate  and 
desorption  rate  have  to  coincide  in  order  to  prevent  the  system 
to  become  over  specified.  Hence,  the  overall  process  is  limited  by 
the  sorption  rate.  The  molar  water  fraction  fv  in  above  equations  is 
taken  to  be  constant  at  0.31.  This  is  justifiable  considering  that  the 
specific  surface  area  is  varied  over  one  order  of  magnitude  in  this 
study. 
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The  equilibrium  concentration  CWiequii  was  determined  by 
Zawodzinski  et  al.  [34]: 

nmem  nmem 

CW,equil  =  WMw  X  =  WMw[(0'°43  +  17  810 

-39.85a2  +36.0a3)sg  +  16.8s;]  (12) 

where  a  is  the  activity  of  water  vapor  adjacent  to  the  membrane 
and  sg  and  S/  are  the  gas  phase  and  liquid  phase  saturation  (volume 
fraction),  respectively.  To  summarize,  the  only  transport  equation 
solved  in  computational  Domain  IV  is  Eq.  (8).  It  is  important  to 
realize  that  the  source  terms  are  non-zero  only  inside  the  catalyst 
layers.  In  the  membrane  the  right  hand  side  is  zero  and  dissolved 
water  transport  reduces  the  diffusive  transport  without  a  source 
term. 

Of  note  is  also  that  there  is  an  enthalpy  term  associated  with  the 
sorption/desorption  of  water  into  and  out  of  the  electrolyte  phase, 
although  this  term  has  often  been  neglected  in  computational  fuel 
cell  models.  Sorption  of  water  vapor  into  the  electrolyte  is  similar  to 
condensation,  while  desorption  is  similar  to  evaporation,  thus  if  the 
overall  drag  is  positive,  i.e.  from  anode  to  cathode,  we  have  a  heat¬ 
ing  term  at  the  anode  and  a  cooling  term  at  the  cathode.  The  heat  of 
absorption  and  desorption  has  been  experimentally  determined  to 
be  very  close  to  the  heat  of  condensation/evaporation  of  water  [35]. 
This  energy  source/sink  term  is  accounted  for  in  Domain  I.  One  fun¬ 
damental  difference  between  the  water  diffusing  through  the  MEA 
in  Domain  IV  and  the  convection  term  calculated  in  Domain  V  is  the 
enthalpy  of  phase  change  which  is  not  accounted  for  when  the  liq¬ 
uid  water  is  merely  pressed  through  the  membrane  by  the  liquid 
phase  pressure  gradient  between  cathode  and  anode.  This  model 
is  hence  also  fundamentally  different  from  the  approach  used  by 
Eikerling  et  al.  [36]  who  considered  EOD  and  hydraulic  permeation 
inside  the  membrane  due  to  capillary  pressure  forces. 

3.4.  Domain  V 

This  domain  accounts  for  the  hydraulic  permeation  of  water 
through  the  membrane.  The  membrane  is  considered  a  porous 
medium  and  the  liquid  phase  pressure  difference  between  cath¬ 
ode  and  anode  side  drives  the  water  through  the  water-filled  pores 
in  the  membrane  [37].  Consequently,  the  only  equation  of  interest 
is  the  three-dimensional  momentum  equation  for  the  liquid  phase: 

V  ■  e(p,V,  ®  U,  -  M,(VU,  +  (VU/)r))  =  e  (  U,  -  Vp,)  (13) 


The  partial  molar  volume  of  the  dry  membrane  is  defined  as: 
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The  partial  molar  volume  of  water  is  simply: 


Vw  = 


Mh20 

Pi 


(17) 


where  MH2o  is  the  molecular  weight  of  water  (18kgkmol-1)  and 
Pi  is  the  water  density  (974  kg  m-3).  The  local  water  content  X  is 
taken  from  the  corresponding  control  volume  in  Domain  IV. 

The  pressure  boundary  conditions  applied  here  result  in  a  cal¬ 
culated  flux,  and  this  flux  in  turn  is  implemented  as  source/sink 
term  for  liquid  water  in  the  CL  cells  adjacent  to  the  membrane  in 
computational  Domain  I.  Note  that  in  all  cases  investigated  here 
the  convective  water  flux  was  at  least  two  orders  of  magnitude 
lower  than  the  fluxes  calculated  in  Domain  IV  and  could  have  been 
neglected.  However,  this  part  of  the  model  was  kept  for  the  sake  of 
completeness  and  because  in  the  future  it  will  allow  us  to  investi¬ 
gate  the  water  flux  when  a  pressure  difference  between  anode  and 
cathode  side  is  applied. 


4.  Model  results 

The  model  was  applied  to  study  the  water  balance  for  a  fuel 
cell  under  specified  operating  conditions.  The  operating  tempera¬ 
ture  was  80  °C  in  every  case  investigated  and  the  outlet  pressure 
was  fixed  at  1.5  atm.  In  stoichiometric  flow  ratio  was  set  to  10 
in  order  to  keep  the  study  at  a  quasi  two-dimensional  level  and 
simplify  postprocessing.  It  should  be  stressed  that  the  model  is 
three-dimensional  in  nature  without  any  numerical  restrictions. 
For  this  particular  study  the  predicted  water  balance  was  inves¬ 
tigated  for  five  different  cases.  Case  1  employed  an  inlet  relative 
humidity  of  100%  at  both  anode  and  cathode  side,  Case  2  an  inlet 
RH  of  50%/50%  and  Case  3  was  conducted  for  entirely  dry  (0%  RH) 
inlet  conditions.  Cases  4  and  5  were  conducted  for  relative  inlet 
humidities  of  0%/100%  C/A  and  100%/0%  C/A,  respectively. 

4.1.  General  remarks 

This  study  was  conducted  in  order  to  shed  additional  light  into 
the  transfer  mechanisms  of  water  across  the  membrane  for  a  fuel 
cell  during  operation.  A  key  output  parameter  is  the  effective  drag 
of  water  across  the  membrane,  typically  defined  as  (e.g.  [20]): 


As  was  discussed  in  our  previous  work  [7]  for  small  bulk  veloci¬ 
ties  U /  the  left  hand  side  of  the  equation  becomes  negligible  and  the 
right  hand  side  reduces  to  Darcy’s  law.  The  computational  domain 
consists  of  the  membrane  only,  where  again  there  is  a  one-to-one 
correspondence  between  these  computational  cells  and  the  mem¬ 
brane  computational  cells  in  Domain  I.  Boundary  conditions  are 
prescribed  in  form  of  two  pressure  boundaries,  where  the  specified 
pressures  correspond  to  the  liquid  phase  pressure  in  the  membrane 
adjacent  CL  cells  from  Domain  I.  The  permeability  of  the  membrane 
is  a  function  of  the  membrane  water  content  calculated  in  Domain 
IV: 


kreiKsat  — 


(14) 


where  I<sat  is  the  liquid  saturated  permeability  (=1.8  x  10-18m2) 
[38].  The  first  part  is  the  relative  permeability  and  it  depends  on 
the  local  wetting  status  of  the  membrane  X  according  to  [37]: 
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where  nw  an  in  denotes  the  incoming  molar  water  stream  at  the 
anode  and  nw_an>0ut  the  outcoming  molar  stream  of  water  at  the 
anode  in  [mols-1].  The  difference  is  of  course  the  net  amount  of 
water  transferred  across  the  membrane.  This  amount  is  normalized 
by  a  half  amount  of  water  produced,  where/'  denotes  the  nominal 
current  density  of  the  operating  cell  in  [A cm-2],  A  is  its  geomet¬ 
rical  area  in  [m2]  and  F  is  the  Faraday  constant  (96,485  C  mol-1). 
A  positive  effective  drag  coefficient  means  that  we  have  net  water 
transport  from  anode  to  cathode,  and  a  negative  effective  drag  coef¬ 
ficient  means  that  we  have  back  transport  from  cathode  to  anode. 
In  the  latter  case  the  back  diffusion  outweighs  the  EOD.  In  our  study 
we  generally  employed  an  EOD  coefficient  of  unity,  so  that  a  result¬ 
ing  effective  drag  coefficient  of  unity  means  that  there  is  no  net  back 
diffusion.  In  this  study  we  can  neglect  convective  transport  of  water 
due  to  the  low  pressure  difference  employed  so  that  we  reduce  the 
overall  water  transport  to  two  phenomena:  EOD  and  back  diffu¬ 
sion.  As  for  diffusion  it  can  be  said  that  the  diffusivity  depends  on 
the  water  content  of  the  membrane:  the  larger  the  water  content  X 
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Fig.  2.  (a)  Predicted  net  drag  coefficient  for  Case  1  as  function  of  the  specific  elec¬ 
trolyte  surface  area  a  [46].  (b)  Predicted  average  membrane  water  content  X  for  Case 
1  as  function  of  the  specific  electrolyte  surface  area  a  [46]. 


the  larger  is  the  membrane  water  diffusion  coefficient.  Hence  a  bet¬ 
ter  humidified  membrane  should  yield  a  lower  (or  more  negative) 
net  drag  coefficient  compared  to  a  dryer  membrane  under  other¬ 
wise  similar  conditions.  We  assumed  a  constant  EOD  of  unity  under 
all  conditions;  the  results  will  indicate  that  A,  >3  for  all  conditions 
investigated  which  verifies  this  simplification  [2]. 

During  this  study  it  soon  emerged  that  the  specific  surface  area 
a  of  the  electrolyte  in  the  CL  plays  a  predominant  role  on  the  overall 
water  transport.  As  this  value  is  very  difficult  to  determine  in  prac¬ 
tice  and  it  can  be  easily  modified  simply  for  varying  the  amount 
of  electrolyte  added  to  the  catalyst  ink  adjacent  to  the  membrane, 
we  conducted  this  study  for  three  different  values  of  a  which  are 
5.0  x  104  m2  nrr3  (“low  a”),  1.0  x  105  m2  nrr3  (“medium  a”),  and 
5.0  x  105  m2  nr3  (“high  a”),  respectively. 

4.2.  Cases  without  MPL 

This  study  contains  cases  with  and  without  an  MPL.  In  the  first 
section  the  MPL  was  left  out.  The  overall  thickness  of  the  GDL  was 
the  same  whether  it  contained  an  MPL  or  not. 

4.2.1.  Casel 

In  the  fully  humidified  case  there  should  be  only  net  water  trans¬ 
port  from  anode  to  cathode  because  there  is  a  near-zero  RH  gradient 
across  the  membrane  which  means  that  the  driving  force  for  back 
diffusion  is  low.  Fig.  2a  and  b  shows  the  resulting  net  drag  coef¬ 
ficient  and  the  membrane  water  content  as  function  of  current 
density  for  the  different  specified  values  of  the  specific  electrolyte 
surface  area  a  in  the  CL.  The  drag  coefficient  is  fairly  independent 
of  the  current  density  which  is  in  agreement  with  experiments  [20] 
but  it  varies  strongly  with  a.  The  membrane  water  content  is  high¬ 
est  for  the  highest  specified  value  of  a,  and  this  should  result  in  a 
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Fig.  3.  (a)  Predicted  MEA  water  content  X  for  Case  1  at  a  current  density  of  0.8  A  cm-2 
and  for  a  low  specific  electrolyte  surface  area  a  =  5.0  x  104  m2  rrr3.  The  cathode  CL 
is  the  lower  part  and  the  anode  CL  the  upper  [46].  (b)  Predicted  MEA  water  content 
X  for  Case  1  at  a  current  density  of  0.8  A  cm-2  and  for  a  high  specific  electrolyte 
surface  area  a  =  5.0  x  105  m2  rrr3.  The  cathode  CL  is  the  lower  part  and  the  anode  CL 
the  upper  [46]. 


superior  cell  performance.  Hence  for  these  operating  conditions  a 
water  uptake  layer  should  be  beneficial. 

Fig.  3a  and  b  shows  the  predicted  membrane  water  content  as 
calculated  in  Domain  IV  at  a  current  density  of  0.8  A  cm-2,  and  they 
demonstrate  an  additional  reason  why  the  performance  should  be 
highest  for  the  case  with  high  a.  The  water  content  inside  the  MEA 
exhibits  a  strong  variation,  and  it  is  interesting  to  note  that  the 
anode  side  water  content  is  predicted  higher  than  the  cathode  side 
while  it  is  common  presumption  that  the  anode  side  of  the  mem¬ 
brane  is  prone  to  dehydration.  While  the  anode  water  content  is  at 
maximum  (A  ^  13)  the  cathode  CL  water  content  is  a  strong  func¬ 
tion  of  a.  It  will  be  seen  below  that  the  fundamental  reason  for 
this  behavior  is  the  variation  of  the  net  drag  coefficient  with  a.  In 
all  cases  the  net  transport  is  from  anode  to  cathode  (see  Fig.  2). 
Even  though  we  encounter  significant  liquid  water  saturation  at 
the  cathode  CL  of  above  30%  (see  Fig.  4)  the  electrolyte  phase  is 
not  fully  hydrated,  which  is  in  contrast  to  the  frequently  assumed 
equilibrium  between  electrolyte  phase  and  gas  (liquid)  phase.  In  all 
of  the  contour  plots  the  results  have  been  mirrored  so  that  the  full 
land  area  is  shown  in  the  middle  (-0.5  mm  to  0.5  mm)  of  the  plot 
and  two  half-channel  areas  at  the  sides. 

Finally,  Fig.  4  shows  that  the  predicted  liquid  water  content  (cal¬ 
culated  in  Domain  I)  at  the  cathode  side  is  only  weakly  dependent 
on  both  current  density  and  a,  but  it  is  very  interesting  to  note  that 
the  anode  side  water  content  shows  a  strong  dependence:  while 
the  anode  always  remains  in  the  single  phase  at  low  current  den¬ 
sities  and  for  a  low  a  it  becomes  multi-phase  at  a  high  current 
density  and  high  a.  Fig.  4  shows  that  under  the  latter  conditions 
there  is  a  significant  amount  of  water  predicted  at  the  anode  side 
under  land.  This  water  content  of  around  20%  corresponds  to  the 
irreducible  saturation,  meaning  that  it  is  not  transported  by  capil¬ 
lary  action.  The  interface  between  anode  GDL  and  channel  remains 
dry,  no  liquid  water  enters  the  anode  flow  channel.  This  water  will 
probably  not  impact  the  cell  performance  but  it  will  be  detected  in 
Neutron  Radiography  experiments. 

For  the  sake  of  completeness  the  relative  humidity  distribution 
inside  the  cell  (Fig.  5a  and  b)  and  the  local  evaporation  rates  (Fig.  6a 
and  b)  are  shown  as  well.  Both  are  calculated  at  a  current  density 
at  0.8  A  cm-2  for  the  low  a  case  and  the  high  a  case,  respectively. 
The  relative  humidity  changes  drastically  at  the  anode  side  when 
we  apply  a  different  a.  While  we  observe  a  dryer  gas  phase  at  the 
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Fig.  4.  (a)  Predicted  liquid  saturation  at  the  cathode  (lower)  and  anode  (upper)  half 
cell  for  Case  1  at  a  current  density  of  0.8  A  cm-2  and  a  =  5.0  x  104  m2  rrr3  [46].  (b) 
Predicted  liquid  saturation  at  the  cathode  (lower)  and  anode  (upper)  half  cell  for 
Case  1  at  a  current  density  of  0.8  A  cm-2  and  a  =  5.0  x  105  m2  rrr3  [46]. 
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Fig.  5.  (a)  Predicted  relative  humidity  distribution  at  the  cathode  (lower)  and  anode 
(upper)  half  cell  for  Case  1  at  a  current  density  of  0.8  A  cm-2  and  a  =  5.0  x  1 04  m2  irr 3 
[46].  (b)  Predicted  relative  humidity  distribution  at  the  cathode  (lower)  and  anode 
(upper)  half  cell  for  Case  1  at  a  current  density  of  0.8  A  cm-2  and  a  =  5.0  x  1 05  m2  irr 3 
[46]. 
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Fig.  6.  (a)  Predicted  local  rate  of  phase  change  at  the  cathode  (lower)  and  anode 
(upper)  half  cell  for  Case  1  at  a  current  density  of  0.8  A  cm-2  and  a  =  5.0  x  1 04  m2  itt3  . 
(b)  Predicted  local  rate  of  phase  change  at  the  cathode  (lower)  and  anode  (upper) 
half  cell  for  Case  1  at  a  current  density  of  0.8  A  cm-2  and  a  =  5.0  x  105  m2  rrr3. 


anode  due  to  the  positive  drag  coefficient  in  the  case  of  low  a  the 
anode  side  remains  fully  humidified  for  high  a. 

Fig.  6  shows  that  the  local  evaporation  rate  inside  the  cathode 
porous  media  becomes  larger  and  the  condensation  zone  becomes 
smaller  with  increasing  a.  The  condensation  zone  near  the  CL  is  due 
to  the  oxygen  consumption,  as  was  described  by  Berning  and  Djilali 
[39].  To  summarize,  it  was  observed  above  that  an  increasing  a  in 
this  case  means  that  the  net  drag  becomes  smaller,  i.e.  lower  water 
transport  to  the  cathode.  The  consequence  of  this  is  that  removing 
liquid  water  from  the  cathode  will  be  less  problematic. 

4.2.2.  Case  2 

In  this  case  the  inlet  gases  were  partly  humidified  at  both  half 
cells.  Due  to  the  water  production  at  the  cathode  this  results  in 
the  build-up  of  a  significant  concentration  gradient  between  anode 
and  cathode.  This  concentration  gradient  is  further  increased  by 
the  EOD.  Hence  compared  to  the  fully  humidified  case  there  is  a 
stronger  driving  force  for  back  diffusion,  and  this  can  even  cause 
the  overall  water  transport  from  cathode  to  anode,  i.e.  a  negative 
net  drag  coefficient. 

Fig.  7a  shows  that  in  this  case  the  net  drag  coefficient  can  either 
be  positive  or  negative,  depending  on  the  specific  surface  area  of  the 
electrolyte  in  the  CL.  A  high  value  of  a  increases  the  back-diffusion. 
The  net  drag  generally  increases  with  current  density  meaning  that 
the  effect  of  the  EOD  becomes  relatively  stronger.  Fig.  7b  shows  that 
the  average  water  content  inside  the  membrane  is  quite  similar  in 
all  cases  at  values  around  A  =  7-8.  A  side  effect  of  using  a  water 
uptake  layer  is  that  we  can  avoid  the  liquid  regime  at  the  cathode 
CL  altogether,  which  leads  to  an  increased  limiting  current  density, 
as  shown  in  Fig.  7c.  Here  we  see  from  the  average  oxygen  concen¬ 
tration  that  for  low  a  the  cathode  CL  is  in  the  two-phase  region 
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Fig.  7.  (a)  Predicted  net  drag  coefficient  for  Case  2  as  function  of  the  specific  electrolyte  surface  area  a  [46].  (b)  Predicted  average  membrane  water  content  k  for  Case  2  as 
function  of  the  specific  electrolyte  surface  area  a  [46].  (c)  Predicted  average  oxygen  concentration  in  the  CL  for  Case  2  as  function  of  the  specific  electrolyte  surface  area  a 
[46].  (d)  Predicted  average  relative  humidity  in  the  CL  for  Case  2  as  function  of  the  specific  electrolyte  surface  area  a. 


for  all  current  densities  and  in  the  single  phase  region  for  high  a, 
while  it  undergoes  a  transition  between  0.6  A  cm-2  and  0.8  A  cm-2 
for  medium  a.  Fig.  7d  summarizes  the  average  relative  humidity  in 
anode  and  cathode  CL  as  function  of  a.  A  water  uptake  layer  signif¬ 
icantly  lowers  the  humidification  gap  between  anode  and  cathode 
and  can  prevent  or  reduce  cathode  flooding,  as  was  observed  above. 
Moreover,  it  may  help  to  “recycle”  water  when  anode  and  cathode 
operate  in  counter-flow  and  the  gases  enter  at  a  fairly  low  RH  as 
was  proposed  by  Biichi  and  Srinivasan  [40]. 

From  Fig.  7a  it  can  also  be  observed  that  at  a  current  density  of 
1 .0  A  cm-2  the  net  drag  is  positive  for  low  a  and  negative  for  high  a. 
Fig.  8a  and  b  shows  the  detailed  calculated  distribution  of  the  elec¬ 
trolyte  water  content  X  for  these  cases.  In  case  of  a  positive  net  drag 
(Fig.  8a)  the  concentration  gradient  of  water  inside  the  membrane 
points  from  anode  to  cathode,  i.e.  the  cathode  side  is  dryer  than 
the  anode  side,  while  for  a  negative  net  drag  the  concentration  gra¬ 
dient  points  from  cathode  to  anode.  This  is  a  direct  consequence 
of  Eq.  (8),  which  describes  the  water  transport  in  the  electrolyte 
membrane  as  pure  diffusion.  This  implies  that  it  depends  entirely 
on  the  net  drag  of  water  through  the  membrane,  which  side  will 
become  dryer.  While  it  is  common  perception  that  the  membrane 
dries  out  on  the  anode  side  (because  the  gas  phase  here  is  dryer), 
these  results  indicate  that  it  is  just  as  likely  that  the  cathode  side 
membrane  may  be  the  dry  side.  Although  this  finding  is  a  direct 
consequence  of  Eq.  (8),  it  does  not  agree  with  experiments  con¬ 
ducted  by  Biichi  and  Scherer  who  found  that  it  is  always  the  anode 
side  that  becomes  dry  [41  ]. 

The  cathode  side  of  the  fuel  cell  has  much  larger  mass  trans¬ 
port  losses  than  the  anode  side,  and  as  a  consequence  the  local 
current  density  moves  away  from  the  CL/membrane  interface  with 
increasing  current  density  [29].  A  drier  electrolyte  phase  inside  the 
cathode  CL  means  that  there  will  be  increased  protonic  loss  here 


and  this  should  be  prevented  in  order  to  optimize  the  cell  perfor¬ 
mance.  Consequently,  operating  conditions  that  cause  a  high  net 
drag  from  anode  to  cathode  should  be  avoided.  A  water  uptake 
layer  may  also  prevent  such  condition. 

Finally,  Fig.  9a  shows  the  predicted  liquid  saturation  inside  the 
cell  at  a  current  density  of  1 .0  A  cm-2  for  a  low  value  of  a.  This  water 
can  be  completely  removed  by  increasing  a  and  hence  reversing 
the  net  drag  coefficient.  There  is  no  liquid  water  under  the  same 
conditions  for  high  a  (not  shown),  and  the  predicted  RH  distribution 
stays  just  below  100%,  as  shown  in  Fig.  9b.  As  was  observed  above, 
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Fig.  8.  (a)  Predicted  membrane  water  content  k  for  Case  2  at  a  current  density  of 
1.0  A  cm-2  for  the  low  specific  electrolyte  surface  area  a  [46].  (b)  Predicted  mem¬ 
brane  water  content  A  for  Case  2  at  a  current  density  of  1.0  A  cm-2  for  the  high 
specific  electrolyte  surface  area  a  [46]. 
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Fig.  9.  (a)  Predicted  liquid  saturation  at  the  cathode  (lower)  and  anode  (upper)  half 
cell  for  Case  2  at  a  current  density  of  1.0  A  cm-2  for  a  =  5.0x  104m2  irr3  [46].  (b) 
Predicted  RH  distribution  at  the  cathode  (lower)  and  anode  (upper)  half  cell  for 
Case  2  at  a  current  density  of  1.0  A  cm-2  for  a  =  5.0  x  105  m2  m-3  [46]. 


this  significantly  increases  the  predicted  limiting  current  density 
from  around  1 .0  A  cm-2  at  a  lower  a  to  over  1 .2  A  cm-2  at  high  a. 

4.2.3.  Case  3 

In  order  to  simplify  system  cost  and  complexity  it  is  desirable 
to  operate  a  low  temperature  PEMFC  without  external  humidifica¬ 
tion.  However,  it  is  commonly  believed  that  this  leads  to  membrane 
dehydration  and  consequently  to  a  low  cell  performance.  In  early 
works  Biichi  and  Srinivasan  experimentally  investigated  the  possi¬ 
bility  to  operate  the  fuel  cell  without  external  humidification  [40]. 
Later,  Biichi  and  Scherer  built  fuel  cells  with  varying  membrane 
thickness  and  found  that  the  limiting  current  density  and  over¬ 
all  membrane  resistance  is  strongly  dependent  on  the  membrane 
thickness  and  hence  on  the  ability  of  the  membrane  to  transport 
water  to  the  anode  via  back  diffusion  [41  ]. 

Case  3  was  carried  out  for  the  case  of  dry  incoming  reactant 
gases.  We  assume  that  such  a  case  can  only  function  with  a  nega¬ 
tive  net  drag  coefficient,  as  otherwise  there  will  be  no  water  at  the 
anode  side  to  sustain  the  EOD.  Hence,  the  limiting  current  density 
is  reached  when  the  net  drag  approaches  zero.  Fig.  10  shows  the 
predicted  development  of  the  net  drag  coefficient.  The  drag  coeffi¬ 
cient  is  in  the  negative  regime  and  approaches  zero  with  increasing 
current  density  for  all  cases.  In  case  of  a  high  a  the  net  drag  remains 
negative  throughout  the  entire  current  density  range  investigated. 

If  we  work  from  the  hypothesis  that  the  limiting  current  density 
is  reached  when  the  net  drag  becomes  zero,  then  for  the  low  a  case 
the  overall  cell  current  is  limited  for  the  water  transport  through 
the  membrane  at  0.5  A  cm-2  and  for  the  medium  a  case  at  around 
1.0  A  cm-2.  Applying  a  water  uptake  layer  may  mean  that  the  lim¬ 
iting  current  density  is  not  determined  by  the  ability  of  the  cell 
to  shed  water  to  the  anode.  The  predicted  membrane  water  con¬ 
tent  is  shown  in  Fig.  10b.  Despite  the  fact  that  the  gases  enter  the 
cell  completely  dry  the  membrane  water  level  is  around  A  «  5  and 
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Fig.  10.  (a)  Predicted  net  drag  coefficient  for  Case  3  as  function  of  the  specific  elec¬ 
trolyte  surface  area  a  [46].  (b)  Predicted  average  membrane  water  content  A  for 
Case  3  as  function  of  the  specific  electrolyte  surface  area  a  [46].  (c)  Predicted  aver¬ 
age  relative  humidity  inside  the  CL  for  Case  3  as  function  of  the  specific  electrolyte 
surface  area  a. 


varies  little  with  current  density.  There  is  only  a  very  small  varia¬ 
tion  in  A  predicted  as  function  of  the  specific  surface  area.  While  a 
partly  dehydrated  membrane  will  lead  to  a  worse  cell  performance 
due  to  high  ohmic  losses  inside  the  membrane  it  may  be  still  worth 
considering  for  some  applications  due  to  the  significant  reduction 
in  system  complexity  and  cost  (no  humidifier  required).  Because 
the  net  drag  coefficient  is  always  negative  in  this  case  this  means 
that  the  electrolyte  phase  in  the  cathode  CL  is  better  humidified 
than  in  the  anode  CL.  Moreover  there  are  no  liquid  water  problems 
expected  at  this  set-up,  but  this  may  change  for  realistic  stoichio¬ 
metric  flow  ratios.  Fig.  11a  and  b  shows  the  predicted  membrane 
water  content  for  the  medium  a  case  and  the  high  a  case  at  a  current 
density  of  1.0  A  cm-2,  where  the  cathode  CL  is  significantly  better 
humidified  for  the  high  a  case  due  to  the  high  back  diffusion  while 
there  is  again  a  near-zero  gradient  for  the  medium  a  case. 
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Fig.  11.  (a)  Predicted  average  membrane  water  content  k  for  Case  3  at  a  current  den¬ 
sity  of  1 .0  A  cm-2  and  a  medium  specific  electrolyte  surface  area  a  [46].  (b)  Predicted 
average  membrane  water  content  k  for  Case  3  at  a  current  density  of  1 .0  A  cm-2  and 
a  high  specific  electrolyte  surface  area  a  [46]. 


We  can  conclude  from  this  that  a  high  specific  surface  area  of 
electrolyte  may  yield  acceptable  currents  when  operating  the  cell 
under  dry  inlet  conditions.  It  is  advisable  to  conduct  experiments 
under  these  operating  conditions  using  different  MEAs  with  vary¬ 
ing  electrolyte  loading  in  the  CL  ink  in  order  to  verify  these  findings. 

4.2.4.  Case  4 

The  fourth  case  was  performed  with  a  fully  humidified  anode 
( 1 00%  RH)  and  completely  dry  air  entering  the  cathode  (0%  RH).  This 
was  also  done  in  order  to  demonstrate  the  numerical  robustness 
of  our  model.  In  addition,  there  are  publications  of  experimental 
works  that  make  this  case  suitable  for  verification. 

In  general  the  results  of  this  case  are  quite  comparable  with  Case 
1,  where  both  anode  and  cathode  were  fed  with  fully  humidified 
gases.  Fig.  12a  shows  that  the  resulting  drag  coefficient  was  pre¬ 
dicted  positive  for  all  value  of  a  investigated,  but  it  decreased  with 
current  density,  which  means  that  relatively  speaking  the  effect  of 
back  diffusion  becomes  stronger  with  current  density.  This  in  turn 
is  in  part  due  to  the  better  humidified  membrane  (see  Fig.  12b), 
which  leads  to  a  higher  diffusion  coefficient,  and  in  part  due  to  an 
increase  in  RH  at  the  cathode  CL  due  to  higher  water  production, 
which  increases  the  driving  force  for  diffusion.  A  detailed  analysis  of 
the  results  revealed  that  all  cases  were  in  the  single  phase  regime  at 
both  anode  and  cathode  and  the  anode  RH  inside  the  CL  was  always 
higher  than  the  cathode  RH.  Note  again  that  the  net  water  transport 
can  be  positive  due  to  EOD  even  when  the  RH  inside  the  cathode 
CL  is  higher  than  the  anode  RH.  The  predicted  membrane  water 
content  increases  with  current  density  from  values  around  X  «  9  to 
around  X  « 12.  As  the  net  water  transport  was  positive  throughout 
this  case  the  cathode  side  membrane  was  predicted  dryer  than  the 
anode  side  membrane. 

Experimental  verification  of  some  of  these  findings  can  be  found 
in  publications  by  Watanabe  et  al.  [42]  who  placed  Platinum  wires 
inside  a  Nation  membrane  to  measure  the  resistance  during  fuel 
cell  operation.  They  found  that  the  membrane  resistance  decreased 
with  increasing  current  density  under  similar  operating  conditions 
(fully  humidified  anode,  dry  cathode).  We  can  confirm  this  trend 
in  our  simulations.  On  the  other  hand  Biichi  and  Scherer  [43]  mea¬ 
sured  the  membrane  resistance  of  Nation  117  by  pulse  response, 
which  however  included  the  electrolyte  resistance  in  the  CL’s  as 
well.  They  used  similar  experimental  conditions  and  found  that  the 
overall  resistance  increased  with  current  density  which  appeared 
to  contradict  the  results  of  Watanabe  et  al.  [42].  To  clarify  this  we 
conducted  simulations  using  a  thicker  (180  pan)  membrane  that 
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Fig.  12.  (a)  Predicted  net  drag  coefficient  for  Case  4  as  function  of  the  specific  elec¬ 
trolyte  surface  area  a  [46].  (b)  Predicted  average  membrane  water  content  k  for  Case 
4  as  function  of  the  specific  electrolyte  surface  area  a  [46]. 


resembled  Nation  117.  It  should  be  kept  in  mind  that  under  such 
operating  conditions  with  a  fully  humidified  anode  and  a  dry  cath¬ 
ode  inlet  gas  the  effective  drag  will  invariably  be  positive,  which 
means  that  the  cathode  CL  will  be  dryer  than  the  anode  CL.  Fig.  13 
shows  the  calculated  membrane  water  content  as  well  as  the  aver¬ 
age  water  content  in  the  cathode  catalyst  layer  (CCL).  While  the 
predicted  membrane  water  content  increases  with  current  density 
as  was  also  measured  by  Watanabe  et  al.  [42]  the  CCL  electrolyte 
water  content  decreases  with  low  and  medium  a.  Along  with  the 
fact  that  the  protons  have  to  migrate  deeper  into  the  catalyst  layer 
with  increasing  current  density  (due  to  mass  transport  limitations 
on  the  oxygen  side)  this  may  lead  to  a  non-linear  increase  in  pro¬ 
tonic  resistance  as  observed  by  Biichi  and  Scherer.  Hence  this  model 


Fig.  13.  Predicted  average  membrane  and  cathode  CL  water  content  k  for  a  Nation 
117  membrane  and  inlet  RH  of  100%/0%  (A/C)  as  function  of  the  specific  electrolyte 
surface  area  a. 
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Fig.  14.  (a)  Predicted  net  drag  coefficient  for  Case  5  as  function  of  the  specific  elec¬ 
trolyte  surface  area  a  [46].  (b)  Predicted  average  membrane  water  content  k  for 
Case  5  as  function  of  the  specific  electrolyte  surface  area  a  [46].  (c)  Predicted  aver¬ 
age  oxygen  concentration  in  the  cathode  CL  for  Case  5  as  function  of  the  specific 
electrolyte  surface  area  a  [46]. 


has  helped  to  explain  this  apparent  contradiction  in  these  two 
experiments. 

4.2.5.  Case  5 

In  Case  5  the  cathode  was  fully  humidified  and  the  hydrogen 
entered  completely  dry.  Again  the  hypothesis  was  that  this  cell 
could  only  function  as  long  as  the  net  drag  remained  negative.  This 
case  is  somewhat  comparable  to  Case  3,  where  completely  dry  gases 
entered  on  both  sides.  As  could  be  expected,  the  resulting  net  drag 
coefficient  is  far  in  the  negative  regime  (see  Fig.  14a).  The  limiting 
current  density  for  the  low  a  case  is  around  0.6-0.7Acm-2  while 
for  the  medium  a  and  high  a  case  the  maximum  current  density  is 
not  limited  by  the  capability  to  shed  water  to  the  anode.  Fig.  14b 
shows  that  the  average  membrane  water  content  X  is  yet  again 
fairly  stable  in  all  cases  at  around  7-9,  and  Fig.  14c  suggests  from 
the  oxygen  concentrations  in  the  CL  that  up  to  a  current  density  of 
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Fig.  15.  (a)  Predicted  liquid  saturation  at  the  cathode  (lower)  and  anode  (upper)  half 
cell  for  Case  5  at  a  current  density  of  0.8  A  cm-2  for  a  =  1.0  x  105  m2  m-3  (medium  a) 
[46].  (b)  Predicted  liquid  saturation  at  the  cathode  (lower)  and  anode  (upper)  half 
cell  for  Case  5  at  a  current  density  of  0.8  A  cm-2  for  a  =  5.0  x  105  m2  rrr3  (high  a). 


0.4  A  cm-2  the  cathode  CL  is  in  the  single  phase  regime  for  medium 
o',  while  for  high  a  it  is  in  the  single  phase  regime  up  to  a  current 
of  0.6  A  cm-2. 

Fig.  15a  and  b  shows  the  liquid  saturation  for  two  chosen  cases, 
one  at  medium  a  and  one  for  high  a.  In  the  high  a  case  the  cathode 
flooding  is  significantly  less  pronounced,  and  it  can  also  be  seen  that 
our  code  is  capable  of  predicting  both  a  “dry-to-wet”  and  “wet-to- 
dry”  transition.  The  distribution  of  the  membrane  water  content 
for  these  cases  is  shown  in  Fig.  16a  and  b.  Because  the  net  drag  in 
the  high  a  case  was  much  more  negative  (from  cathode  to  anode) 
the  resulting  gradient  in  water  content  across  the  membrane  is  very 
high,  leading  to  a  well  humidified  cathode  CL.  In  the  medium  a  case 
the  cathode  CL  was  completely  in  the  wet  regime,  as  can  be  seen  in 
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Fig.  16.  (a)  Predicted  membrane  water  content  k  for  Case  5  at  a  current  density  of 
0.8  A  cm-2  and  a  medium  specific  electrolyte  surface  area  a.  (b)  Predicted  membrane 
water  content  k  for  Case  5  at  a  current  density  of  0.8  A  cm-2  and  a  high  specific 
electrolyte  surface  area  a. 


6316 


T.  Berning  et  al.  /  Journal  of  Power  Sources  196(2011)  6305-6317 


LIQ  SAT  [-] 


o 


5.0E-04  - 

-1.0E-03  -5.0E-04  0.0E+00  5.0E-04  1.0E-03  1.5E-03 

Cell  width  [m] 

Fig.  17.  Predicted  liquid  saturation  at  the  cathode  (lower)  and  anode  (upper)  half 
cell  for  Case  1  at  a  current  density  of  0.8  A  cm-2  for  a  =  5.0  x  105  m2  m-3  when  an 
MPL  was  included  at  anode  and  cathode  side. 


Fig.  1 8.  Predicted  net  drag  coefficient  for  Case  1  as  function  of  the  specific  electrolyte 
surface  area  a  when  an  MPL  was  included. 


Fig.  1 5a,  but  the  resulting  water  content  in  the  electrolyte  was  only 
around  9-10,  hence  far  from  its  equilibrium  value  (see  Fig.  16a). 

4.3.  Cases  with  MPL 

When  an  MPL  was  included  it  was  30  fxm  thick.  As  listed  in 
Table  1  the  MPL  was  modeled  having  a  lower  permeability,  a  higher 
contact  angle  by  comparison  to  the  GDL  and  CL  and  a  lower  irre¬ 
ducible  saturation  while  the  porosity  was  assumed  to  be  the  same 
as  for  the  GDL.  The  tortuosity  was  also  increased  for  the  MPL.  Note 
that  in  a  CFD  model  the  porosity  predominantly  plays  a  role  in 
transient  simulations;  it  cancels  out  of  the  equations  in  steady 
state  simulations  (except  for  the  species  diffusion  coefficients).  The 
above  changes  to  the  material  properties  mean  that  the  MPL  is 
more  hydrophobic  and  less  permeable  than  the  GDL,  and  it  provides 
a  larger  diffusion  resistance  due  to  the  higher  tortuosity.  These 
changed  material  parameters  have  also  been  applied  in  a  similar 
manner  by  Pasaogullari  and  Wang  [44]. 

While  all  cases  as  above  were  conducted  in  the  same  detail  with 
MPL  we  can  summarize  the  results  by  stating  that  the  predicted 
impact  of  the  MPL  on  the  water  balance  and  membrane  water  con¬ 
tent  was  very  small  for  all  cases.  In  most  cases  there  appeared  to 
be  no  difference  at  all  in  the  predicted  results  whether  an  MPL  was 
included  or  not. 

Fig.  17  shows  the  predicted  liquid  water  saturation  for  Case  1 
at  a  current  density  of  0.8  A  cm-2  for  a  high  a  when  an  MPL  is 
included.  This  compares  to  Fig.  4b,  where  no  MPL  was  included. 
We  can  clearly  see  the  predicted  saturation  jumps  between  CL  and 
MPL  and  between  MPL  and  GDL  [7].  As  was  discussed  by  Nam  and 
Kaviany  [27]  these  jumps  arise  out  of  the  condition  that  the  liquid 
phase  pressure  (capillary  pressure)  is  continuous  over  the  inter¬ 
faces,  but  the  saturation  versus  capillary  pressure  functions  vary 
for  each  layer.  Pasaogullari  and  Wang  had  similar  findings  [44]. 
However,  when  looking  at  the  predicted  net  water  drag  for  this 
case  including  an  MPL  we  can  see  that  there  is  almost  no  effect 
predicted.  Fig.  18  shows  very  similar  results  to  Fig.  2a. 

The  strongest  effect  of  the  MPL  was  found  in  Case  5,  where  the 
cathode  gas  enters  fully  humidified  and  the  anode  gas  enters  com¬ 
pletely  dry.  While  the  predicted  water  balance  is  very  similar  with 
and  without  MPL  (compare  Figs.  19a  and  14a)  addition  of  the  MPL 
leads  to  a  higher  limiting  current  density  in  the  case  of  high  a  (com¬ 
pare  Figs.  19b  and  14c).  This  is  owing  to  the  fact  that  significantly 
less  liquid  water  is  predicted  at  the  cathode  side  at  elevated  current 
densities  (compare  Figs.  20  and  15b). 

Overall  it  must  be  stated  that  in  general  the  addition  of  an 
MPL  in  our  model  by  modifying  material  parameters  such  as 
apparent  contact  angle  (increased  hydrophobicity),  permeability 
(increased  pressure  drop),  irreducible  saturation  (smaller  fraction 
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Fig.  19.  (a)  Predicted  net  drag  coefficient  for  Case  5  as  function  of  the  specific  elec¬ 
trolyte  surface  area  a  when  an  MPL  was  included,  (b)  Predicted  average  oxygen 
concentration  in  the  cathode  CL  for  Case  5  as  function  of  the  specific  electrolyte 
surface  area  a  when  an  MPL  was  included. 
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Fig.  20.  Predicted  liquid  saturation  at  the  cathode  (lower)  and  anode  (upper)  half 
cell  for  Case  5  at  a  current  density  of  0.8  A  cm-2  for  o' =  5.0  x  105  m2  m-3  (high  a) 
when  an  MPL  was  included. 
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of  hydrophilic  pores)  and  tortuosity  (increased  diffusion  resistance) 
did  not  noticeably  affect  the  water  balance  or  membrane  water  con¬ 
tent.  This  is  in  accord  with  experimental  work  conducted  by  Karan 
et  al.  [22]. 

5.  Conclusions 

We  presented  a  CFD  model  of  a  PEMFC  that  has  novel  fea¬ 
tures  and  utilizes  the  multifluid  approach.  The  computational 
domain  was  split  up  into  five  subdomains  which  greatly  improved 
convergence  behavior  at  the  cost  of  having  a  slightly  higher  num¬ 
ber  of  computational  cells.  A  study  was  carried  out  under  quasi 
two-dimensional  operating  conditions  (short  active  area,  high  stoi¬ 
chiometric  flow  ratios)  in  order  to  gain  fundamental  understanding 
of  the  transport  processes  involved.  The  following  conclusions  can 
be  drawn  from  this  study: 

•  The  overall  water  balance  can  depend  strongly  on  the  spe¬ 
cific  surface  area  of  electrolyte  in  the  catalyst  layers  because 
for  thin  membranes  the  net  water  transport  might  be  limited 
by  sorption/desorption  kinetics.  This  can  potentially  be  verified 
experimentally  by  varying  the  electrolyte  loading  close  to  the 
CL/membrane  interface. 

•  The  micro-porous  layer  was  not  found  to  have  a  notable  impact 
on  the  water  content  or  water  transport  across  the  membrane. 
In  accordance  with  previous  studies  it  was  found  that  the  water 
saturation  level  inside  the  MPL  is  lower  than  inside  the  GDL,  but 
this  does  not  impact  the  predicted  saturation  level  inside  the  CL, 
nor  does  it  affect  the  net  water  balance.  More  work  is  needed  on 
fundamentally  understanding  the  role  of  the  MPL  in  addition  to  its 
beneficial  effect  due  to  decreased  electronic  resistance  between 
the  GDL  substrate  and  the  CL  [45]. 

•  The  net  drag  coefficient  alone  defines  which  side  of  the  electrolyte 
membrane  becomes  dehydrated.  This  was  concluded  from  fun¬ 
damental  considerations  by  looking  at  the  transport  equation  for 
water  inside  the  membrane  which  reduces  to  a  diffusion  equation 
with  no  source  terms.  Consequently,  a  high  positive  drag  (from 
anode  to  cathode)  means  that  there  is  a  large  concentration  drop 
of  water  inside  the  membrane  from  anode  to  cathode,  i.e.  the 
cathode  becomes  dry.  This  state  should  be  avoided  as  then  the 
electrolyte  phase  in  the  cathode  CL  is  predicted  dry,  too,  which 
means  that  we  encounter  high  protonic  losses.  At  high  current 
densities  the  local  current  generation  in  the  CL  moves  away  from 
the  membrane  due  to  mass  transport  limitations  which  further 
amplifies  the  detrimental  effect  of  a  dry  electrolyte  phase. 

•  If  we  start  from  the  hypothesis  that  the  cell  can  only  operate  when 
the  anode  side  is  prevented  from  drying  out  we  may  be  able  to 
operate  the  cell  under  dry  anode  inlet  conditions  by  adjusting  the 
specific  surface  area  of  the  electrolyte  in  the  CL. 

•  Even  when  operating  the  cell  under  completely  dry  conditions 
the  membrane  water  content  was  predicted  to  be  4-5,  hence  far 
from  being  completely  dry. 

•  On  the  other  hand  the  membrane  water  content  was  not  always 
at  its  maximum  when  both  gases  enter  the  cell  fully  humidified 
(Case  1 ).  A  high  net  water  drag  means  that  the  cathode  side  of  the 
membrane  dries  out. 

•  Anode  liquid  water  was  found  only  at  high  current  densities  and 
for  high  values  of  the  specific  surface  area  of  the  electrolyte  in  the 
CL.  This  water  was  located  under  land  only  and  was  not  predicted 
to  reach  the  GDL/channel  interface. 


•  Flydraulic  permeation  of  water  through  the  PEM  was  found  to  be 
at  least  two  orders  of  magnitude  smaller  than  diffusion/EOD. 

Most  of  these  findings  can  be  verified  by  conducting  suitable 
experiments,  which  will  be  the  focus  of  our  future  work. 
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